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REFRESHED SHOCKS 
and AFTERGLOW LONGEVITY in GRB 
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OS : ABSTRACT 



QJ , We consider fireball models where the ejecta have a range of bulk Lorentz factors 

so that the inner (lower T) parts may carry most of the mass, or even most of the 
OO I energy. The outer shock and contact discontinuity decelerate as the fireball sweeps 

up external matter. This deceleration allows slower ejecta to catch up, replenishing 
and reenergizing the reverse shock and boosting the momentum in the blast wave. In 
^ ' consequence, the energy available to power the afterglow may substantially exceed that 

'^ ■ of the burst itself. Such models allow a wide range of possibilities for the afterglow 

^vq . evolution, even in the case of spherically symmetric expansion. 



Subject headings: gamma-rays: bursts 



1. Introduction 

Afterglows from Gamma-Ray Bursts (GRB) have been discovered in at least seven objects 
at X-rays wavelengths. Two of these also yielded optical afterglows and one was followed in up 
to the radio band (e.g. papers at the 4th Huntsville GRB Symposium, Meegan et al, 1997). 
The first afterglow discovered at both X-ray and optical wavelengths, GRB 970228 (Costa et 
al, 1997) strongly favored a cosmological origin, and this was strengthened even further by the 
discovery in the second object, GRB 970508, of several systems of absorption lines yielding a limit 
0.835 ^ z ^ 2.3 (Metzger et al, 1997). The power law fiux decay predicted from the simplest 
cosmological fireball afterglow model (Meszaros & Rees, 1997) turned out to be in good general 
agreement with the observed behavior (e.g. Tavani, 1997| ; |Waxman, 1997a ; Vietri, 199"7b|jWijers^ 



Rees fc Meszaros , 1997;Reichart, 1997]). The simple model is also successful in explaining 



the rise and decay of the light curves ([Wijers, Rees &: Meszaros , 1997|), as well as the overall 



radio behavior ( Goodman, J., 1997 ; [Waxman, et al, 1997c ). The ensuing massive observational 



campaigns on these sources have, in the meantime, provided such excellent data that it is now 
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feasible to investigate more subtle effects. Among these are, e.g. possible oscillations superposed 
on the power law decay (e.g. Galama, et al., 1997), the appearance of optical and X-ray bumps 



after ~ 1.5 days on the optical and X-ray light curves ( piro, et al, 1997 ), the evidence for an 



initial slow optical decay before the onset of the power law rise and decay (Pedersen, et al, 1997), 
etc. Most and perhaps all of these features may be understood within the context of anisotropic 
outflows, or of isotropic outflows in inhomogeneous media ( pVIeszaros , Rees fc Wijers 1997 ). It 



remains, however, an interesting question to investigate how slow can a power law decay behavior 
be in a simple spherically symmetric model, and what fraction of the total energy budget can 
be associated with the late afterglow, as opposed to the initial GRB. In most earlier discussions, 
the ejecta have been treated as a uniform sphere or shell, where most of the energy is carried by 
material with a well-deflned (and generally high) Lorentz factor. But it is possible - perhaps, 
indeed, more physically realistic - that the ejecta consist, in effect, of many concentric shells 
moving at different (relativistic) speeds. This could be the outcome, for instance, of relativistic 
shock propagation down a density gradient; or the inner shells might be more heavily loaded with 
baryonic ejecta. A given shell would then expand freely until the contact discontinuity had been 
decelerated, by sweeping up external material, to a Lorentz factor lower than its own. It would 
then crash into a reverse shock, thermalizing its energy and boosting the power of the afterglow. 
In this paper, we illustrate this phenomenon by discussing simple cases. We show that such models 
can produce very slow decay rates in the afterglow; they can accommodate very large (as well as 
small) ratios of afterglow to burst fluences, and can naturally explain a variety of afterglow light 
curves. 



2. Kinematics of Ejecta with a Power-Law Radial Profile 

The 'trigger' for a gamma ray burst may have complex time structure spread over the 
duration (typically 10 seconds) of the intense gamma ray emission itself. However, when we are 
considering mechanisms for the (much more prolonged) afterglow, we can treat this energy release 
as impulsive. It gives rise to relativistically expanding debris which eventually slows down as it 
sweeps up external matter. If the debris were spherical, it would be characterized by the way its 
energy and mass were distributed among shells of different Lorentz factor. The actual T-profile' 
would be the outcome of complex dynamics during the burst itself: shock fronts propagating down 
density gradients, internal shocks, etc. A range of Lorentz factors seems likely - indeed it is not 
obvious whether the energy (or even the mass) would be concentrated towards the high or low 
end of the T distribution. This is our motivation for first exploring a range of power-law' models, 
which can be treated analytically. In simple cases where the expansion is ultrarelativistic (bulk 
Lorentz factor T ^ 1) and all physical variables depend as power laws on radius (or source-frame 
time), the dynamics can be treated in a 1-D approximation, which will be valid either for spherical 
expansion or for jet angles ^ T~^. We consider the evolution at large radii, after the input of 
mass and energy from the central source has ceased. In the source frame, after deceleration by 
the external medium becomes important, the Lorentz factor of the outer shock and the contact 



discontinuity begin to decrease as 

Tc cc r-, (1) 

When the ejecta have, in effect, a single (high) Lorentz factor, we have the well-known models 
where n = 3/2 for adiabatic (weak coupling), or n = 3 for radiative (strong coupling) expansion in 
a homogeneous external medium (for an inhomogeneous medium, see Meszaros , Rees & Wijers, 
1997). In the cases treated here, where the ejecta has a range of Lorentz factors Fj, the expansion 
will still obey a power law, but with different n. At increasing r the contact discontinuity, if it 
expands according to (1), lags behind the light cone by a fractional amount 



oc — / r^^(r) dr oc — / r"^" dr oc r^^(r) . (2) 

r 2r J "" ^ ' 2r J 2(2n + 1) '^ ^ ^ ^ ' 

We assume that the input of energy is variable, but occurs over a finite time which is short 
compared to the deceleration time of the contact discontinuity, so the impulsive approximation is 
valid. Subsequent shells of fireball material endowed with a decreasing Lorentz factor Tf, catch up 
with the contact discontinuity at radii r which satisfy (Ar/r) = [2(2n + l)]~"'^r~^ = (l/2)r7 , or 

r^(r) = (2n+l)i/2r,(r), (3) 

where we have neglected the thickness of the shell of reverse shocked gas relative to r. This is 
shown in Figure 1, where r is distance from the center of the outflow and t is source- frame (not 
detector) time, t = r/c. The impacting material therefore has a Lorentz factor which is {2n + 1)^'^ 
times larger than that of the contact discontinuity. (The latter has been decelerated, whereas the 
former, having moved at constant speed, reaches the same radius with a higher speed, see Figure 
1). The ratio of the Lorentz factors, however, is only a factor ^^ 2, so the reverse shock(s) are 
marginally relativistic at all times. It then makes little difference to the overall dynamics whether 
the shock is radiative or adiabatic. 



3. Effects of a Fj-dependent Mass and Energy Input 

We assume that the primary event leading to the GRB produces a fireball in which the mass 
fraction ejected above a given initial bulk Lorentz factor decreases according to 

M(> F/) oc rj' . (4) 

The mass is dominated by the low Tf shells for s > 0, while the energy will also be dominated 
by low Tf material if s > 1. The mass that has 'caught up' with the contact discontinuity, and 
passed through the reverse shock, before it gets out to radius r is (assuming Ar/r <C 1) 

M(< r) oc r"^ (5) 

The momentum impacting on the contact discontinuity from the reverse shocked gas, per unit 
solid angle and measured in the comoving frame, is 

Pre. OC '-^^ OC r^^ . (6) 

dr 



The impact is always marginally relativistic, Trev ~ 2 ~ constant. The force on the other (forward 
shock) side of the contact discontinuity is 

Pfor OC Pext^l^^r"^ , (7) 

where ^ = for radiative (strong coupling) and A = 1 for adiabatic (weak coupling) expansion of 
the remnant (Meszaros , Rees & Wijers, 1997), and the geometrical factor r^ is included because 
this is per unit solid angle. Here the adiabatic or radiative dependence is important, since the 
outer shock is ultrarelativistic. We note that ^ = and A = \ are extremes of the remnant 
behavior; intermediate values are also expected in realistic cases, e.g., as shown in the numerical 
calculations of Panaitescu & Meszaros , 1997. For pext independent of radius, equating the powers 
of r in equations (^) and ^ gives ns — 1 = — n(l + A) + 2, so the index n = 3/(s + 1 + ^) and 
the radius dependence of the contact discontinuity Lorentz factor in the presence of a continuous 
input of mass and energy is 

Fc ex r"" ex r-3/(^+i+^) . (8) 

The value of this index is shown in Table 1 for various cases. 

We see that, for s = (where the later arriving, slower mass makes only a logarithmic 
difference to the total), we recover as expected the previous law of n = 3/2(3) for the A = 0(1) 
cases corresponding to the adiabatic (radiative) dynamic regimes. However for s = 1, where most 
of the mass is concentrated at lower Fj we get n = 1(3/2) for A = 0(1); and for s = 3/2, where 
also most of the energy is dominated by the low Fj, we get n = 6/7(6/5) for A = 0(1). 



4. Luminosity and Flux Evolution 

For s > 1 the total energy available for the afterglow gradually increases during the 
deceleration, but clearly the power law dependences will only be valid over a finite range of 
times. Keeping these limitations in mind, the luminosity dependence at late times can be 
evaluated. For a contact discontinuity evolution given by equation (^, the energy liberated by 
the shock is E cc FJ~^ oc 7.-3(i-s)/(s+i+A) ^^^ ^.j^^ observer time over which this is released is 
T oc r/Tl oc r('^+^+'^)/('^+^+^), corresponding to r oc t'^s+i+A)/{s+A+7) ^ ^pj^g "kinetic" luminosity or 
flux is then 

where T is detector (observer-frame) time. E.g., for s = 3/2, A = 1 we have Fj^in oc T~^^''^'^ . 

If the radiation comes from the forward shock, and is due to synchrotron radiation of 
shock-accelerated power law electrons in a turbulently generated magnetic field which has built 
up to some fraction of equipartition with the protons, we have a comoving synchrotron peak 
frequency v'^ oc ^^B' oc F|? oc r~^/('^+^+^) oc 7"~9/('5+^+7)^ ^nd the corresponding observer-frame 
value is Vm oc Fcf^ oc 7"-i2/(s+A+7)^ where 7e is the electron random Lorentz factor. Following 



Meszaros , Rees fc Wijers 1997 the comoving synchrotron peak intensity is I^, oc n'^-yei''^ e, 



sy, 



where n^ oc TcUext ~ ^c is the comoving postshock electron density for expansion in a homogeneous 
external medium, and Csy oc (i^wn/^sy)" ^ r'^F^" is the synchrotron radiative efficiency of the 
electrons responsible for radiating at the peak frequency. For a = the electrons are radiatively 
efficient and Cgy = 1, while for a = 1 the electrons are adiabatic (radiatively inefficient) and 
Csy < 1. Thus II^ oc rT^^-i oc T^<^is+i+A)-6a+3]/{s+A+7) _ rpjjg observer frame peak synchrotron 
flux is 

F,„ oc T^Tlll^ oc T-l^a-is+l+A)(2+a)]/{s+A+7) (^^q) 

For a synchrotron intensity spectral slope P above the synchrotron peak, I^ (x u^ , the time decay 
of the flux in a fixed detector frequency band i^d (e-g- optical) is 

Fd oc F,,Jz.b/j.^)^ oc r[(2+«)(«+l+^)-6«+12/3]/{^+^+7) (i^) 

For s = 3/2, A = 1 this is oc T^7-(5a/2)+i2p]/{i9/2) ^^^^ for a = 1, /? = 1/3 (e.g. synchrotron 
radiation below the peak from adiabatic electrons) the flux would increase oc T^^'^^, steeper than 
in our previous models. For a = 1, f3 = —1 (above the peak) the flux is oc T~^^'^^. Similar 
considerations can be made for the reverse shock; however, if Rayleigh- Taylor instabilities occur 
across the contact discontinuity (e.g. Waxman & Piran, 1994), these would tend to equalize the 
magnetic and nonthermal energy densities, leading to a similar time-dependence of the emission. 
Further complications may arise if the electrons responsible for Vm are adiabatic but those 
responsible for z^o are radiative, e.g. equation (9) in [Meszaros , Rees fc Wijers 1997 and related 



discussion (also |Sari, Piran fc Narayan, 1997| ). Note that equation (11) above was calculated 
assuming /3 ^ — 1. For P <> —1 most of the energy does not reside at Um but at some largest 
frequency Ucut to which the slope /3 extends to, and the equation replacing ( |ll|) would depend on 
this cutoff frequency. The cutoff can be computed, given model hypotheses defining the absolute 
values of the magnetic field and Tc- Note also that equations (^,PJ1^, 11) can be easily generalized 
to take into account a dependence on an external medium that depends on radius as next oc r~'^, as 
in Meszaros , Rees &: Wijers 1997, leading to flatter decay laws than in equation (0). In Table 1 



we show the values of the time exponents of the fluxes in a fixed detector band for a homogeneous 
external medium and some representative values of s. A, a and (3. 



5. Discussion 

We have considered here GRB models where the ffieball does not constitute a single 
ultrarelativistic shell (or uniform sphere), but contains material with a range of speeds, larger 
amounts of mass and energy being characterized by lower bulk Lorentz factors. After the central 
energy-generation has ceased the impulsive approximation is valid, and during its evolution the 
later portions of the ejecta catch up with the previously decelerated matter, thus continually 
re-energizing the shock that provides the afterglow of the GRB. We have shown in §^ that, 
depending on the Lorentz factor dependence of the mass and energy injection, even spherically 
symmetric outflows can produce very flat flux time decay laws (equations ||ll|). Even flatter 



decay laws are possible for synchrotron energy spectral indices P ^ —1. A significant departure 
from our previous versions of an isotropic GRB afterglow model is that here the flux decay law 
is not dependent only on the spectral index (or on the density or angular dependence indices, 
c.f. Meszaros , Rees & Wijers, 1997) but also on the index s of the Lorentz factor dependence of 
the mass and energy ejection (equation Q). Even for a homogeneous medium and for isotropic 
outflow, plausible values of s can lead to afterglow luminosities whose time-integrated value can 
either exceed or fall below the initial 7-ray fluence. The former case would represent a situation 
similar to the hypernova model of Paczyiiski, 1997. 

Although we have focused on power-law dependences (cf equation (1) and (4)), our 
considerations clearly show that even a spherical model could generate more complicated light 
curves. For example, suppose that, rather than containing material with a smooth distribution 
of Tf, the fireball could be better modeled by two shells with very different Tf, the slower 
shell carrying more energy. The afterglow luminosity would initially decline as in the standard 
mono-energetic case; but when the blast wave had slowed down enough for the slower shell to 
catch up, there would be a boost to the luminosity before it again resumed a power-law decline. 
Changes in the slope of the flux decay laws (equations |ll[),or late bumps (e.g. Piro et al, 1997), 
could thus be readily incorporated even into spherical models. The "internal engine" associated 
with a catastrophic disruption or collapse event would inevitable lead to non-sphericity in the 
outflow; if the ejecta consisted of discrete overdense blobs, their deceleration could lead to more 
complex time-structure, or apparent shorter-term variability (Galama et al, 1997, Pedersen, et al, 
1997) superposed on the overall power law behavior. 
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Fig. 1. — Schematic space-time diagram in source frame coordinates of a relativistic outflow, with 
a range of Lorentz factors Tf, triggered by an explosion that can be approximated as instantaneous 
(decreasing Tf values lead to world lines further to the left, which are straight before entering 
the shock). The contact discontinuity and the forward shock are being decelerated due to the 
increasing amount of external matter being swept up, so that they lag behind the light cone by 
an increasing amount amount Ar whose increase with r is steeper than linear. This deceleration 
allows slower ejecta (with lower Tf) to catch up and pass through a reverse shock just inside the 
contact discontinuity. Both shocks are here approximated as having negligible width. Energy 
and momentum from progressively slower material is thus continuously re-energizing the shocks, 
resulting in a more gradual decay of the afterglow. 
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Table 1: For strong coupling ("radiative", ^4 = 0) or weak coupling ( "adiabatic" , A = 1) remnant 
evolution in a homogeneous external medium and a mass outflow index s {M{> Tf) oc T7^) the 
bulk Lorentz of the contact discontinuity evolves as Tc oc r~", producing a kinetic flux F^j^ ex t^ 
whose index v is in column 5. For different spectral indices /3 the forward shock produces a flux 
in a fixed detector band Fd oc t"", whose indices w are shown in columns 6 through 8 for radiative 
(a = 0) or adiabatic (a = 1) electrons. Flatter decay laws would be obtained if the external density 
decreases with radius. 



